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Abstract 

The rotating string model is an effective model of mesons, in which the quark and the antiquark 
are linked by a straight string. We previously developed a new framework to include the retardation 
effects in the rotating string model, but the string was still kept straight. We now go a step further 
and show that the retardation effects cause a small deviation of the string from the straight line. We 
first give general arguments constraining the string shape. Then, we find analytical and numerical 
solutions for the string deformation induced by retardation effects. We finally discuss the influence 
of the curved string on the energy spectrum of the model. 
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I. INTRODUCTION 



The retardation effect between two interacting particles is a relativistic phenomenon, due 
to the finiteness of the interaction speed. Light mesons are typical systems in which these 
effects can significantly contribute to the dynamics, since the light quarks can move at a 
speed close to the speed of light. We developed in Ref. [l| a generalization of the rotating 
string model (RSM) in order to take into account the retardation effects in the mesons. 
The RSM is an effective model derived from the QCD Lagrangian, describing a meson by a 
quark and an antiquark linked by a straight string. Both particles are considered as spinless 
because spin interactions are sufficiently small to be added in perturbation. It has been 
shown that the RSM was classically equivalent to the relativistic flux tube model 



This last model, firstly presented in Refs. ja, |7[, yields meson spectra in good agreement 
with the experimental data 5]. Our method to treat the retardation effects relies on the 
hypothesis that the relative time between the quark and the antiquark must have a nonzero 
value. Consequently, in our approach, the evolution parameter of the system is not the 
common proper time of the quark, the antiquark and the string, but the time coordinate of 
the center of mass which plays the role of an "average" time. 

We showed in Ref. jl| that, in the special case where the quark and the antiquark have 
the same mass, the part of the total Hamiltonian containing the retardation terms could be 
treated as a perturbation. This perturbation is a harmonic oscillator in the relative time 
variable, with an effective reduced mass and an effective restoring force both depending on 
eigenstates of the unperturbed Hamiltonian (which is independent of the relative time). The 
fundamental state of this oscillator gives the contribution of the retardation to the masses 
as well as the relative time part of the wave function. So, the relative time wave function is 
a gaussian function centered around zero. This point confirms the validity of the usual non 
retarded model in first approximation. It is worth mentioning that our retardation term 
does not destroy the Regge trajectories, the linear relation between the square mass and the 
spin of the light mesons. Our generalized RSM also allows to reproduce the experimental 
meson spectrum with a good agreement. 

This previous work must be considered as a first trial to compute an estimation of the 
retardation effects in mesons. Some hypothesis were made in order to keep the calculations 
workable. In particular, the straight line ansatz was used to describe the string. Although it 
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simplifies the calculations, it is worth noting that the use of a non vanishing relative time is 
not really compatible with a straight string. We already gave in Ref. |l| a crude estimation 
of the possible bending of the string due to retardation effects, and our result was compatible 
with a small deviation from the straight line. However, this point deserves a further study to 
confirm the validity of our approach. Moreover, these deformations are of intrinsic interest, 
since they can exist independently of the retardation effects. 

Our paper is organized as follows. In Sec. Ull we briefly recall the approach developed 
in Ref. ^, the RSM with a nonzero relative time. Then, assuming that the string can be 
curved, we give arguments constraining its possible shape in Sec. IIIH and we make a rather 
general ansatz for the curved string in Sec. IIVI Using this ansatz, we obtain analytic and 
numerical solutions for the string shape in Sees. |3 and IVII Finally, we sum up our results 
in Sec. IVTll 



II. ROTATING STRING MODEL WITH A NONZERO RELATIVE TIME 

The RSM with a nonzero relative time has been studied in detail in Ref. So, we 
simply recall here the main points of this work. Starting from the QCD Lagrangian and 
neglecting the spin contribution of the quark and the antiquark, the Lagrange function of a 
meson reads [3] ('7 = diag(H ) and h = c = 1) 

C = —miy xl — m2y xl — a / d6 y {ww')'^ — w'^w''^ . (1) 
^ Jo 

The first two terms are the kinetic energy operators of the quark and the antiquark, whose 
current masses are mi and m2- These two particles are attached by a string with a tension a. 
Xi is the coordinate of the quark i and w is the coordinate of the string, w depends on two 
variables defined on the string world sheet: One is spacelike, 6', and the other timelike, r. 
Derivatives are denoted w' = dgw and w = drW. In this picture, r is the timelike evolution 
parameter for the string and the quarks. Let us mention that bold quantities will always 
denote four-vectors. 

Introducing auxiliary fields to get rid of the square roots in the Lagrangian (^J and making 

the straight line ansatz to describe the string, an effective Lagrangian can be derived 

1 
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2{c2 + Ca2)rr + a^r' + {a4 + 2Cci + Cai)r' , (2) 
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where the coefficients Oi , 02, . . . , are functions of the auxihary fields. Their exact expressions 
can be found in Ref. The parameter ( defines the position R of the center of mass 

R = Cxi + {l-C)x2 = {t,R), (3) 
and r is the relative coordinate 

r = Xi- X2 = {(T,r). (4) 

It is worth mentioning that the auxiliary field /ij can be interpreted as the constituent mass 
of the quark whose current mass is rrii. Moreover, the parameter ( appears to be a function 
of the auxiliary fields. It reduces to the expected value ( = 1/2 when the quark and the 
antiquark have the same mass . The straight line ansatz for the string implies that the 
string coordinates are given by 

w = R+{e -C)r. (5) 

Such an ansatz is suggested by lattice QCD calculations, which show that the chromoelectric 
field between the quark and the antiquark appears to be roughly constant on a straight line 
joining the two particles j^. 

The usual approach is to work with the equal time ansatz, i. e. 

a;0 = x° = = r = i. (6) 

Then, we have r = (0,r), R = (t, R ), r = (0,r), and R = (1,-R). This procedure 
considerably simplifies the equations, but neglects the relativistic retardation effects due to 
a possible nonzero value of the relative time a. That is why we made in Ref. [l| a less 
restrictive hypothesis: We identified the temporal coordinate of the center of mass with the 
evolution parameter, t = t, and we allowed a non vanishing relative time a. We have then 

r = {&,n, R={l,n). (7) 

It is then possible to derive from the Lagrangian Q a set of three equations for the RSM 
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(8a) 
(8b) 

(8c) 
(9) 



arcsin?/i -y,^l-yf 

Pr is the radial momentum and is the transverse velocity of the quark i. The first relation 
gives the cancellation of the total momentum in the center of mass frame, while the two 
last ones define respectively the angular momentum and the Hamiltonian of the system. 
Equations ()8a|) and ()8b|) define the two variables yi and ?/2- A further elimination of the 
auxiliary fields /ii and /i2 yields the relativistic flux tube model 

Equations (jHl) are identical to those of the usual RSM (see for example Ref. P]), but a 
perturbation of the Hamiltonian, denoted AH, is now present. It contains the contribution 
of the retardation effects and is given by 

C2 • 

- -I So- - (ci + Cai)a - 

as Zas 



AH 



2a'- 



1 • • 

+ — Sa - (ci + Qai)a - + -(04 + 2Cci + Cax)ci^ 



(10) 



where S is the canonical momentum associated with the relative time a. 

Let us now consider the quantized version of our model: L ^ + 1)? Vr 
[a, S] = —i. The Hamiltonian (jHcj) has then the following structure 



H{a,r) = Ho{r) + AH{a,r). 



(11) 



The relative time only appears in the perturbation, and Hq depends only on the radius r. 
So, we can assume that the total wave function reads 



{r)) = \Aia))0\Rir))0\YeUOA)), 



(12) 



where \R{r)) is a solution of the eigenequation 



Ho{r) \R{r)) = Mo\R{r)) 



(13) 
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Such a problem can be solved for instance by the Lagrange mesh technique As it is 
shown in Ref. the total mass is written 

M = Mo + {A{(y) I ® {R{r) \ AH{r, a) \R{r)) ® \A{a)) 

= Mo + AM. (14) 

The contribution AM is then given by the fundamental state of the eigenequation 

AHia) \A{a)) = AM \A{a)) , (15) 

where 

An{a) = {R{r) \ AH{r, a) \R{r)) . (16) 

In the case mi = m2, the eigenequation ()15|) can be solved, since A7i (cr) is simply a harmonic 
oscillator in the relative time a 

A^(a) = --^ [S^ + {cl - a,a,)a^] . (17) 

A more detailed study of the retardation term AM and of its effect on the meson spectrum 
can be found in Ref. Let us only observe two points from Eqs. fjTHjl and (fTTjl . Firstly, 
the retardation contribution AM is negative. 



y{cl-a,as)/{as)^ (18) 



AM 

2 

and it thus decreases the meson mass. Secondly, the temporal part of the wave function 
reads 

1/4 



A{a)={^ exp(-^a^), (19) 



with 



P = ^J{cl-a^as). (20) 



It is a gaussian function centered around a = 0. This provides an interpretation of the equal 
time ansatz (jH)) as the most probable configuration of the system. 

It is worth noting that the use of a non vanishing relative time is not really compatible 
with the straight line ansatz. This can be seen by the following simple considerations: Let 
us assume that the world sheet of the system in the center of mass frame is a helicoid area 
in the case of exactly circular quark orbits. The shape of the string is then a straight line 
for a slice at constant time and a curve for a slice not at constant time (see Fig. We 



FIG. 1: Left: Intersection between a helicoidal string world sheet and a plane r = tq {x^{tq) = 
x^{tq)), the intersection is a straight line. Right: Same situation, but now Xg(ro) 7^ x^{tq), and 
the intersection is curved. 

already gave in Ref. [l| a crude estimation of the possible bending of the string due to 
retardation effects, and our result was compatible with a small deviation from the straight 
line. The purpose of this paper is to study more carefully the string deformations caused by 
the nonzero relative time in order to confirm the validity of our approach. 



III. CONSTRAINTS ON THE STRING SHAPE 

As the coordinates Xi are independent of the parameter 6, the Lagrangian (Q) can be 
rewritten as 

c = [ dec, (21) 

Jo 

with 

jC = —Till \Jd^ — m2\J~^ — a\J {ww'Y — w'^w''^ . (22) 

Several steps are required to obtain a Hamiltonian from the Lagrangian ()22|1 . Firstly, one 
has to find a particular solution, denoted w* , for the string shape. This is in fact the solution 
of the equations of motion (EOM) of the Nambu-Goto Lagrangian 



Cng = -a^J{ww'f - w'^w'^, (23) 

with the boundary conditions 

w{1,t) = x^{t), w{S),t) = X2{t). (24) 



Secondly, once w* is known, it has to be injected in the total Lagrangian ()22|). Thirdly, one 
has to compute the momenta defined by 

In this picture, we assume that w*{6,t) is in fact of the form w*{6, Xi{T)), because of the 
boundary conditions Finally, with the momenta ()25|). the quantity 

2 

n = J2Pi^^-C (26) 

1=1 

is readily computed, and the total Hamiltonian is thus given by 

1 



n= / den. (27) 

Jo 

In the usual RSM, one uses the straight line ansatz (0) together with^ the equal time 
ansatz (jHI). This provides indeed a correct solution w* of the EOM 11, p. 122]. The 
question we want to answer in this section is: Is it possible to find a more general form for 
the string shape? 

Let us suppose that we have such a solution, w*, which thus respects the conditions ((211) • 
As we already remarked, w* depends on the coordinates ^^(t). So, the momenta ()25|) read 

Ptt, = -mi^= - a— ==— ^, (28) 

with 

N = {ww')w' - w'^w. (29) 

For the sake of simplicity, we have dropped the * on w. The Hamiltonian ()2()j) . computed with 
the momenta (j2HI), is expected to vanish since we work in a manifestly covariant formalism. 
Here, we see that 



H = a 



^jNw dx 

If we want = 0, we must satisfy the condition 

2 

dx 



(30) 



ElS^r- (31) 



The solution for the constraint ()31|) is 



«;^(e,r) = ^Af,(^,r)xr(r), (32) 



1=1 



where the A1^{6,t) are unspecified functions. A simple integration gives 

Y.Ai{e,T)x^{T)-Y^ f dfAi{e,f)x^{f) + ne). (33) 



i=l i=l 

Expressing Eq. (I^Hj) in the center of mass and relative coordinates, given by relations ^ 
and (jH), we obtain 

wP = « + Al)R'' + [(1 - CX - CAl] (34) 
- 1 rfr {« + k'^^R' + [(1 - CX - C^L] r'] + 

if we assume that C = 0. We did not write explicitly the dependences in (6*, r) to clarify the 
notations. 

Formula (jMj) can be simplified by taking the limit r — 0. In this case, we want indeed 
to obtain w = R. This clearly implies that 

< + A^, = 5^, g^ = 0, (35) 

and the whole expression can thus be rewritten by using a single set of functions 

w'' = Rf + {AP - CKV - J drAy. (36) 

The solution ()36|) is the general string shape we search for. Solving the corresponding 
EOM would allow us to determine the functions A^^. However, this problem is very complex, 
and we need to make some simplifications to go further. We want to find a solution which 
corresponds to a string configuration which does not change in time. So we consider the 
functions A^ to be independent of r. We also make the following hypothesis 

K = A^K. (37) 

without summation on p. This implies that we want depending only on the coordinates 
R^^ and r'^. With our assumptions, the string has finally the following form 

w^{e, t) = R^{t) + [A^{e) - C] r^(r), (38) 

without summation on /i. The remaining unknown quantities are the functions A^{d). Let 
us note that the straight string ansatz is only a particular case of formula (IHHj). where 

A^'{e) = e. 



IV. THE CURVED STRING ANSATZ 



In order to further simplify the problem, we will here consider mesons in which the quark 
and the antiquark have the same mass. In this case, ( = 1/2, and only the part of the 
string linking the center of mass and one quark has to be known. The second part can be 
computed from the first one by using symmetry arguments, as it will be discussed in Sec. IVII 
Formally, we thus only have to find a solution for w{p, r) with p G [0, 1], such that 

w{0,t)=R, w{l,T)=rgiT), (39) 

rq(r) being the position of the quark, for instance. As we work in the center of mass frame, 
we can set R = (t , ) . We can then rewrite Eq. ()38p on the following form 

w%p,t) = i{T) + A'{py^{T), (40a) 
w\p,t) = A\p)rl{T). (40b) 

The meson evolves in a plane. We can thus set A'^ equal to zero and use the complex 
coordinates {w^,w,w*) defined by |3] 

w = ^{w^ + iw^). (41) 

The Nambu-Goto Lagrangian is invariant for the reparameterization of the world sheet jlll 
p. 93], that is to say the transformations 

T' = r'ip,T), p' = p'ip,T). (42) 

It allows us to fix for simplicity 

t(r) = r, and A^{p) = p. (43) 
will always be assumed. Let us note that t(r) = r was a 



In the following, the relations 
hypothesis of the model in Ref. [l| 

It can be shown that, when r°(r) = 



w°(p, r) = r, (44a) 
w{p,r) = ^^pexp[z^r], (44b) 
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is a solution of the EOM of the Nambu-Goto, or equivalently, of the Polyakov Lagrangian 
(see for example Ref. Q]). In our coordinates, these EOM read 



0, (45a) 



da 

with h""^ the inverse matrix of hab, given by 

hab = daw'^dbW^, 

= -daw^dbw'^ + 2Re{daW*dbw) . (45b) 

The indices (a, 6) label the parameters (t, p), and hab is the induced metric on the string 
world sheet. 

Solution describes a straight string with a constant angular speed, and it corresponds 
to the following choices in Eqs. (j^ : r°(r) = 0, A'^{p) = p, r^i^r) = rqi^r) cosujt, and 
r^(r) = rq(r) sin cut. 

In our previous work on the retardation effects, we used the following ansatz Q 

w°(p, r) = T + pag{T), (46a) 

T ( T ) 

w{p,t) = pexp(zco'r), (46b) 

v2 

that is to say r°(r) = CTq^r), A°(p) = A^(p) = p, rg(r) = rg(r) coscur, and r^(r) = 
rg(r) sin tut. However, the string defined by Eqs. is not a solution of the EOM. In- 
deed, as o"g(r) is assumed to be small (we showed that the retardation can be treated as a 
perturbation), we can write the EOM at the first order in aq. After some calculations, one 
can check from relations ()45p that ag must satisfy the following equation 

ryag - 2rgrqp'ag + [{ry - r^fg + 2rf) - 2] a, = 0. (47) 

As ag is only a function of r, Eq. (jTfj) must be satisfied for every value of p. In particular, 
when p = 0, we observe that = is the only possible solution. This confirms that a 
straight string is not compatible with a nonzero relative time. 

Finally, it appears quite natural that the string could be curved because of the addition 
of two effects: The rotation of the meson and the finiteness of the gluon speed. A curved 
string can be described by the ansatz 

w°(p,r) = r + p(Tg(r), (48a) 
w{p, r) = ^ [p + z/(p)] exp [i{ujr + 0(r)] , (48b) 
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where the spatial deformation /(p) has been introduced as a counterpart to the relative 
time crg(r). Moreover, the eventuality of a non constant rotation speed is taken into account 
through the angular acceleration 4>{t). If ag = 0, our ansatz reduces to the one of Ref. jl^ . 
where the string deformations due to angular acceleration are studied. Equations (jlKjl clearly 
describe a curved string, as it can be seen by rewriting it when r = with 0(0) = 0. Then 
we have simply oc p and w"^ oc /(p). 



V. ANALYTIC SOLUTION 



Finding an exact expression of /(p) which satisfies the EOM of the Nambu-Goto La- 
grangian is a very complex problem, out of the scope of this paper. The use of approxima- 
tions appears necessary in order to deal with workable equations. We have given in Ref. 
some arguments to show that the deformation of the string should be small. Assuming 
that point, we will linearize the EOM in /(p) and its derivatives. The angular accelera- 
tion will also be considered as small, as it is done in Ref. Following the hypothesis of 
this last reference, we will only keep the linear terms in and its derivatives. Terms like 
(pf, (fxTg, (j)dg, (pTq, . . . will thus also be neglected. Moreover, we consider that r^, ~ 0. 
Consequently, our solution will only be valid in the case of small radial excitations. 

After a tedious algebra, we find that the EOM are satisfied by the curved string ()48|l 

if 

(1 _ A,p - A^p' + Asp') dlfip) + (Si - B2p) (pd.fip) - fip)) = C + Dp- Ep\ (49) 
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with 



T 

Ai = (Tq— — 3dq, (50a) 
rq 

A2 = u^l + + 2a,a,^ - 3a/, (50b) 

r 

A3 = Uj'^aqrqVq + (Tg— - C^V^^r^ - dqi-l - (Tqdq^— + dq^ , (50c) 

• 2 

B, = uj^l + 3a;VJ + 2af-\ - 2a,a,^, (50d) 

fq fq 

■'^ • 2 . 



= 2cuV,r,f, + 2a,^ + u'dqT^ - ScuV.a^^ - 2(T,r^^ + 2a'^dq^ + 2(7,^^^^, (50e) 

Tq Tq 
C = 2aqUJ, (50f) 



D = lujoqdq - 2ujaf-^ + r% (50g) 
E = -cjV^rJ + uo^al - 2ujaqrl + 2cuaJ^ + 2ujdqrqrq - 2ujaldq^, (50h) 

'T'q 

and the boundary conditions 

/(0) = /(l)=0. (51) 

Conditions (j^Tj) are in fact equivalent to the initial boundary conditions (jH^ . 

Before performing a numerical resolution of the differential equation ()49|) . we can find an 
approximate analytic solution by developing /(p) in powers of p, 

OD 

/(p) = 5^a„p". (52) 

n=0 

Keeping in this series only the terms which satisfy Eqs. and ()51|) at the second order 
in p, one can find that 

/(p) = /c(p) + /d(p) + /e(p), (53) 

with 

fcip) = -|p(l-p), (54a) 
fnip) = -^-^^^p{l-p\ (54b) 
fE{p) = ^[E-A,{D + CA,) - C{A^ - B,/2)] p (1 - p'). (54c) 

The solution at the lowest order is fc{p)- K is a parabola whose maximum is reached 
in p = 1/2. The next functions, foip) and fsip), shift slightly the maximum at a value 
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p < 1/2. A graphical representation of the solution ()53|) is given in the next section. It is 
worth mentioning that in the case of a vanishing angular momentum, u; = 0, the solution 
is trivially /(p) = 0. Even when the retardation is included, the string is straight when the 
angular momentum is zero. 

We can check that if cTq = 0, our solution reduces to 

/(P)U=0 = -^P(1-P')^ (55) 



which is precisely the result of Ref. 12] • However, as we are only interested in the retardation 



effects, we will take = in the latter. 



VI. NUMERICAL SOLUTION 



Evaluation of the coefficients 



The previous section gave us a qualitative idea of the string shape. But, to complete our 
analysis, we need an estimation of the magnitude of the deformation. To do this, we have to 
compute numerically the different coefficients not in the classical framework we used 
up to now, but in a quantized model. The coefficients (jSUjl should then be seen as operators 
whose average value has to be computed. What should be done rigorously is to rewrite our 
RSM Hamiltonian with a curved string solution of Eq. ()49p and compute the eigenstates 
of this model to average the operators (jSIl), in an analog way of what is done in Ref. jl^j ]. 
But here, we only want to have a first estimation of the string deformation. That is why, as 
we considered that the deformation was small, it appears reasonable to compute the mean 
values of the coefficients (jSUj) with the states of our particular RSM introduced in Ref. 
Let us remark that, if we want to be consistent with the notations for r and a in Sec. m we 
have to make the following substitutions 

r,^r/2, a,^a/2, (56) 

the angular speed u being not modified. Our procedure will be to replace each positive 



definite operator X appearing in the definitions (jHn|) by a/ {X'^), a quantity easy to compute 
with our numerical method |lO]. This simple computation, neglecting symmetrization prob- 
lems, will only give us a rough estimation of the different coefficients, but it is sufficient for 
our purpose. 
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We already mentioned that our numerical method, the Lagrange mesh method, ensures 
us to know the radial wave function \R{r)). The average values (p^) = {p"^ + i{i + l)/r^) 
and (r^) are easy to compute Moreover, it is shown in Ref. [14] that the angular speed 
can be approximated by 



- - (57) 

(r2) (^(^y/p + m^^ + a (r) /g) ' 

where i is the orbital angular momentum. The last spatial term we need is (f^), which is 
the radial part of (r^). It can be computed, from the Lagrangian that j5[ 



(p^) = {{asr - c^rr) . (58) 
Using the fact that (02) ~ [l|, we have 

(r^) - A m - (cl) (r^)] • (59) 
(as) 

We turn now our attention to the terms involving the relative time. With the temporal 
wave function given by Eq. (fTTH) . it is easily computed that 

yp)-^. (60) 



(cT^) can be computed in analogy with (r^). From the Lagrangian (j2I), it can be shown 
that Q 

{^') = {{c,a - a;&f) . (61) 



Since 



we have 



:S^> = f , (62) 



.2\ P 



2 {asf 



1 



(4) 



(63) 



We are now able to compute, at least in first approximation, the coefficients ()5()|1 . 
B. Results 

The more interesting case is the meson formed of two massless quarks. The deformation 
of the string is indeed expected to be maximal in this case since the relativistic effects are 
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TABLE I: Average values of quantities involved in the coefficients (jSOp . for different states. For 



simplicity, the quantities of the form \J {X"^) appearing in the first column are denoted as X. We 

recall that = r/2 and dg = a /2. 

(n+l)L IP IF 2P 

rq (GeV^i) 2.526 3.376 3.411 

(Tq (GeV-^) 1.093 1.072 1.122 

uj (GeV) 0.086 0.087 0.036 

Vq 0.459 0.303 0.592 

a„ 0.342 0.257 0.252 



the most important. Of course, we will choose £ 7^ to observe a nonzero deformation. 
With m = and the standard value a = 0.2 GeV^, we can compute the needed average 
values thanks to formulas ()57|) to ()63p. They are given in Table H] for three different states. 

A numerical integration of Eq. ()49|1 to obtain the solution /(p) with the boundary con- 
ditions (jKT|) can be performed by using the values of Table HI As mentioned in Sec. IIV| one 
half of the string is simply given by the couples {p, /(p)}. Once this solution is known, we 
can compute every couple |^^, /(^)| between the quark and the antiquark. Indeed, a simple 
symmetry argument allows us to write 

r ~ 1 f {(1 + P)/2, /(p)} [1/2, 11, 

^ ^ [{(l-p)/2, -/(p)} ^^G [0,1/2]. 

Relations (jGljl simply define a central symmetry with respect to the center of mass. The 
numerical solution of Eq. pHj) is plotted in Fig. |21 for the three states considered in Table HI 

We found that the deformation is maximal for the IF state, and then decreases when 
the quantum numbers increase. This is a consequence of our approach, since we noticed 
in Ref. |l| that the retardation effects are less large when the dynamical quark mass 
(a/p^ _|_ j77,2^ increases. In every case, the deformation f{9) is small with respect to one. 
This is an a posteriori validation of our choice to consider small deformations only. 

In Fig. ini we compare the numerical solution for the IP state, the lowest state in which 
the deformation occurs, with three analytic approximations of this solution, given by rela- 
tions (fSHj) and It appears that fc is an upper bound of the deformation, and that 
fc + fo is sufficient to correctly approximate the numerical solution. 
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FIG. 2: Numerical solution of Eq. (|49j) . giving the string shape between the quark and the antiquark 
for the states considered in Tabled /(^) = corresponds to the straight line ansatz. 




FIG. 3: Comparison between the numerical solution of Eq. H49|) for the IP state and the analytical 
formula 

Since the string brings an energetic contribution to the meson which is proportional to 
its length, it is interesting to evaluate the ratio between the lengths of both the curved and 
the straight strings. It is given by 



1 



1, 



(65) 
(66) 



For a small deformation, as it is the case here, we can obtain an upper bound for AL/L by 
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using fc instead of /. To the second order in C, we find that 

^ < ^, (67) 

L - 24' ^ ^ 

with C = 2ijqUj. It is maximal in the IF state. We finally obtain the upper bound 



< 2 X IQ-l (68) 

The length of the string is only modified by some tenths of percent, because of the bending 
induced by retardation effects. As the typical mass scale for the mesons is the 1-2 GeV, 
the correction due to curved string is around 1-4 MeV, as it is observed in Ref. The 
contribution of the bending of the string to the mass spectrum seems thus very small. This 
is also small compared with the retardation contribution, which can be around 100 MeV for 
massless quarks. 



VII. CONCLUSION 



We developed in this paper the idea that, if the retardation effects are included in the 
rotating string model, the string linking the quark and the antiquark cannot remain a straight 
line. We found a relation constraining its shape, and obtained a general form for the string. 
The straight line ansatz, which is valid in the equal time approximation, appears to be 
only a particular case of our general form. When the retardation is included, we showed 
that a bending of the string must be taken into account, and we proposed a general ansatz 
defining a curved string. With this ansatz, and assuming a priori that the deformation 
of the string is small, we derived a differential equation giving its shape. Analytical and 
numerical solutions of this equation can be found. Roughly, the deformed string has a 
parabolic shape, with a small amplitude. The amplitude of the deformation decreases when 
the quantum numbers increase; but it is zero for a vanishing angular momentum. We finally 
argued that the contribution of the bending of the string to the mass spectrum is around 
1-4 MeV. The typical meson mass scale is around 1-2 GeV, and the retardation contribution 
computed in a previous work with a straight string is always smaller than 100 MeV. So, the 
effects of the curvature are very weak, although rather interesting from a theoretical point 
of view. However, if one is mainly interested in the computation of a mass spectrum, we can 
conclude from our present work that the rotating string with nonzero relative time can give 
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satisfactory results, event with a straight string approximation. We leave the analysis of 
the rotating string model equations with nonzero relative time and curved string for future 
work. 

Acknowledgments 

C. S. (FNRS Research Associate) and F. B. (FNRS Research Fellow) thank the FNRS 
for financial support. V. M. (IISN Scientific Research Worker) thanks the IISN for financial 
support. 



[1] F. Buisseret and C. Semay, Phys. Rev. D 72, 114004 (2005) |hep-ph/0505T68] . 

[2] A. Yu. Dubin, A. B. Kaidalov, and Yu. A. Simonov, Phys. Atom. Nucl. 56, 1745 (1993); Yad. 

Fiz. 56, 213 (1993) |hep-ph/93H'344l . 
[3] E. L. Gubankova and A. Yu. Dubin, Phys. Lett. B 334, 180 (1994) hep-ph/9408278| . 
[4] C. Semay B. Silvestre-Brac, and I. M. Narodetskii, Phys. Rev. D 69, 014003 (2004) 

|hep-ph/0309256|. 

[5] F. Buisseret and C. Semay, Phys. Rev. D 70, 077501 (2004) |h ep-ph/0406216] . 

[6] D. LaCourse and M. G. Olsson, Phys. Rev. D 39, 2751 (1989). 

[7] M. G. Olsson and S. Veseh, Phys. Rev. D 51, 3578 (1995). 

[8] C. Semay and B. Silvestre-Brac, Phys. Rev. D 52, 6553 (1995). 

[9] Y. Koma, E. M. Ilgenfritz, T. Suzuki, and H. Toki, Phys. Rev. D 64, 014015 (2001) 
|hep-ph/0Q11165| . 

[10] F. Buisseret and C. Semay, Phys. Rev. E 71, 026705 (2005) hep-ph/0409033. . 

[11] B. Zwiebach, A first course in string theory, Cambridge University Press, 2004. 

[12] T. J. Ahen, M. G. Olsson and S. Veseh, Phys. Rev. D 59, 094011 (1999) |hep-ph / 98 10363]. 

[13] T. J. Ahen, M. G. Olsson and S. Veseh, Phys. Rev. D 60, 074026 (1999) |hep-ph /9903222]. 

[14] F. Buisseret and C. Semay, Phys. Rev. D 71, 034019 (2005) |hep-ph/0412361| . 



